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ABSTRACT 

We present a scenario (SupraNova) for the formation of GRBs occurring 
when a supramassive neutron star (SMNS) loses so much angular momentum 
that centrifugal support against self-gravity becomes impossible, and the 
star implodes to a black hole. This may be the baryon-cleanest environment 
proposed so far, because the SN explosion in which the SMNS formed swept 
the medium surrounding the remnant, and the quickly spinning remnant loses 
energy through magnetic dipole radiation at a rate exceeding its Eddington 
luminosity by some four orders of magnitude. The implosion is adiabatic 
because neutrinos have short mean free paths, and silent, given the prompt 
collapse of the polar caps. However, a mass M; 0.1 Mq in the equatorial belt 
can easily reach centrifugal equilibrium. The mechanism of energy extraction is 
via the conversion of the Poynting flux (due to the large-scale magnetic field 
locked into the minitorus) into a magnetized relativistic wind. Occasionally this 
model will produce quickly decaying, or non-detectable afterglows. 


Subject headings: gamma rays: bursts - stars: neutron - black holes - relativity: 
general - instabilities 


1. Introduction 

A successful model for the progenitors of gamma ray bursts (GRBs) must make contact 
with the fireball model (Rees and Meszaros 1992), which so successfully predicted the 
discovery of afterglows in the radio (Paczyhski and Rhoads 1993), in the optical (Meszaros 
and Rees 1997a, Katz 1994), in the X-ray band (Vietri 1997), and explained radio flares 
(Goodman 1997). While some features of afterglows remain at present unexplained, yet 
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there can be little doubt that the overall energetics, time-scales, power-law behaviour, and 
appearance of longer wavelength radiation at later times are all nicely accounted for. 

The crucial point of the hreball model is that it channels the large explosion energy 
(> 10^^ erg, barring beaming) in a locale with small baryon contamination 10“*^ Mq, 
Rees and Meszaros 1992). From this point of view, hypernova explosions (Woosley 1993, 
Paczyhski 1998) are problematic: massive stars are well known to have powerful winds 
which may carry away upward of 10 Mq, thusly baryon-contaminating the environment in 
which the GRB is eventually to go off. Obviously, a very thin beam may be a way out of 
this conundrum, but this begins to look unlikely, given the currently available observational 
evidence (Grindlay 1998). Also scenarios involving accretion-induced collapse of white 
dwarfs (Usov 1992, Blackman, Yi, Field 1996, Yi and Blackman 1997) are troublesome, both 
because of the large amounts of baryons lying about, and because of the long timescales 
required to accrete sizeable amounts of mass. 

Given the considerable specihc angular momenta involved, a plausibly baryon-clean 
model is the merger of neutron star/neutron star, or neutron star/black holes binaries 
(Narayan, Paczyhski and Piran 1992) which is expected to take place outside the 
environment where the binary formed (Bloom, Sigurdsson and Pols 1998). However, at 
least some bursts (though not all, Wijers and Galama 1998) seem to arise within dense 
environments suggesting star forming regions: such is the case of GRB 970111 (Feroci et ai, 
1998), GRB 970828 (Murakami et ai, 1997, Groot et ai, 1998a) and GRB 980326 (Groot 
et ai, 1998b). Since it is possible that GRBs arise from several distinct sources (Meszaros, 
Rees and Wijers 1998b), we propose here a different model, where the initial explosion is 
triggered by the gravitational implosion to a black hole of a supramassive neutron star: we 
call this model a SupraNova. 

In the following, we shall first discuss some properties of equilibrium supramassive 
neutron stars, then the implosion to a back hole and the extraction of the energy that leads 
to the GRB. Reasons why this model remains baryon-clean will be discussed in Section 3, 
and in the last section we will discuss some observational consequences. 


2. On supramassive neutron stars and their implosion 

All known equations of state for matter at nuclear densities allow the construction of 
equilibrium models of supramassive neutron stars (SMNSs), i.e. rotating sequences with 
constant baryon number but different angular momenta; no member of this sequence may 
be non-rotating because they have masses at infinity larger than the largest mass for static 
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models (for recent work on the topic see Cook, Shapiro and Teukolsky 1994, and Salgado, 
Bonazzola, Gourgoulhon and Haensel 1994). Equilibrium sequences for models with baryon 
number exceeding a critical Be above which no static equilibria are possible, start out 
exactly at break-up angular speed; depending upon the equation of state, these models 
have masses in the range M k, 2 — 3.5 Mq, equatorial radii Req ~ 11 — 18 fcm, and angular 
velocities uj ~ 8,000 — 12,000 Though fast rotating, the total angular momentum of 
these models is not large: measured in units of the critical angular momentum R = GM'^/c 
of a black hole of the same mass at inhnity, they have j ~ 0.6 — 0.78. An important feature 
of these models is that, as they lose angular momentum, through most of the models’ 
parameter space they speed up: they contract so as to reduce their moment of inertia, 
and increase their angular velocity. For this reason, angular momentum loss through, for 
instance, magnetic dipole radiation actually speeds up. As the losses mount, the models 
reach a point, having about half of the initial angular momentum, where they become 
secularly unstable to axisymmetric perturbations. The evolution of these models beyond 
this point has not been investigated yet, but the eventual fate of the star is sealed, since it 
cannot connect to a static, stable conhguration: it will implode to a black hole. 

The above results hold for every equation of state, so that we regard them as rather 
robust. Currently, there seems to be no obvious theoretical reason why these stars should 
not form in the collapse of massive stars. A recent claim (Andersson 1998 ) that a new 
type of r-modes leads to such a copious emission of gravitational waves, to make any 
fast-spinning neutron star slow down within a year, would seem to imply that SMNSs, 
which speed up as they lose angular momentum, should collapse to a black hole faster than 
expected on the basis of pure magnetic dipole radiation (Eq. |l|, to be derived shortly). 
However, these computations only concern idealized neutron star models, and neglect 
proper, general relativistic treatment of viscosity, convection, and non-linear effects, all of 
which may slow down, or even damp this instability; they will be neglected henceforth. 

An approximate estimate of the time it takes for this conhguration to collapse is 
given by assuming that the major source of angular momentum loss is through the usual 
magnetic dipole radiation; the star’s rotational energy is Erot = Juj/2, and, since uj is 
roughly constant as the star loses angular momentum passing through different stages of 
equilibrium (Salgado et al, 1994, Cook et al, 1994), Erot ~ ujJ/2; equating this to the 
classical Pacini (1967) formula, one gets J, and the spin-down time to halve J, as required 
by the onset of the instability 


tsd = — = 10 yr — 
J ^ 0.6 
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Since this is much longer than the viscous timescale inside the NS, on which the 
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configuration evolves after it has become secularly unstable toward the hnal collapse, tgd is 
a good estimate of the timescale between the hrst explosion (a normal SN where a SMNS is 
formed) and the implosion to a black hole. 

As the star loses angular momentum, the centrifugal force weakens until support 
against self-gravity becomes impossible: the star implodes. The implosion is likely to 
be strongly non-homologous, the central densest regions forming a black hole hrst, with 
the outermost regions collapsing later. To a hrst approximation, the implosion will be 
‘silent’, with most matter being swallowed down the hole promptly, since the total angular 
momentum of the collapsing NS is subcritical {j ~ 0.6 — 0.78). In particular, it is likely that 
all mass lying close to the spin axis ends up in the hole, a point already made by Meszaros 
and Rees (1997b) and Paczyhski (1998). Furthermore, the neutrinos’ mean free path inside 
the NS is expected to be small, so the collapse will be, most likely, adiabatic. Further 
baryon contamination must be small; for instance, it is difficult to think of any outwardly 
propagating shock capable of ejecting baryons, because of the lack of a hard surface 
and of energetic phenomena capable of releasing considerable counterpressures. Also, for 
reasonable magnetic helds {B < 10^^ G), any magnetic-related phenomena using all the 
energy in the magnetic held could push at most ~ 10“® Mq in centrifugal equilibrium, all 
most likely in the equatorial plane. 

It is likely that this implosion is accompanied by the release of a copious amount of 
energy in gravitational waves; in fact, the known, axisymmetric instability of the equilibrium 
sequence is just secular (Friedman, Ipser and Sorkin 1988); it is possible that this leads to a 
nearby, distinct equilibrium model which becomes dynamically unstable as it loses angular 
momentum; the most likely source of instability is to a bar-formation mode {m = 2) which 
leads to release of gravitational waves (Chandrasekhar 1970). 

Not all mass will be immediately accreted, though: the outermost layers, in fact, have 
centrifugal accelerations only ^ 10% below the local gravitational attraction (Salgado et 
ai, 1994). Thus, even if the pressure gradient were suddenly removed, these layers would 
just contract a bit (barring large energy losses due to shocks caused by the converging 
flow), and halt in centrifugal equilibrium at a radius ~ 10% smaller than the initial one. 
We cannot compute Mi at this stage, so we give only a suggestive argument showing that, 
possibly. Ml ^ 0.1 Mq. Stellar dynamical simulations of self-gravitating galaxies in purely 
rotational equilibrium show that these objects are unstable to bar-forming m = 2 modes, 
a well-known result (Hohl 1971, Ostriker and Peebles 1973, Binney and Tremaine 1987). 
One may see these as energy conserving instabilities, where stars may acquire some energy 
because the gravitational potential in which they move is not constant in time. In such 
systems, typically ~ 1% of the mass (originally, all stars on high angular momentum orbits) 
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reaches orbits outside the main body of the resulting bar, becoming detached (but not 
unbound) from the galaxy. We may expect the total fraction of mass lost in hydrodynamic 
systems to be comparable to this, since here too parcels of fluid on high angular momentum 
orbits interact with the rest of the system mostly through the gravitational potential; in 
fact, according to the models of Salgado et al. (1994), the contribution of the pressure 
gradient to the support of the outer layers of the star is much smaller than the centrifugal 
one. Moreover neutron stars differ from stellar dynamical systems because of losses by 
gravitational waves, and possibly neutrinos, which weaken gravitational binding of the 
outermost layers: it is sufficient that ~ 10% of the whole star mass at inhnity ends 
up as GWs, for the equatorial belt to become rotationally-supported: this corresponds 
to 0.3 Mq radiated away, well within current estimates of GW emission efficiencies in 
non-axisymmetric collapses to black hole (Smarr 1979). Neutrinos not trapped inside the 
collapsing SMNS will also carry away binding mass, or may provide pressure pushing the 
equatorial belt into centrifugal equilibrium. We thus expect that the SMNS may typically 
shed a few percent of its total mass to centrifugal equilibrium. This implies Mi 0.1 Mq. 

The ensuing conhguration is that of a largish black hole (M si 2 — 3 Mq), spinning 
at subcritical speed (j ~ 0.6), surrounded by a thin equatorial belt of matter of mass Mi 
which could not be immediately accreted because its angular momentum, supplemented by 
the factors discussed above, makes it rest in centrifugal equilibrium. The binding energy of 
the left-over mass Mi is 


Eh = 3 X 10 
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This binding energy must be radiated as the leftover accretes onto the black hole, and it 
is this energy release which powers the GRB proper, not the gravitational collapse per se. 
Even the most powerful bursts (Kulkarni et ai, 1998) only require Mi ^ O.W^Mq, where 
6‘^/4 is the beaming fraction. 

The leftover material will be threaded by the same, large-scale magnetic field lines 
that threaded it inside the NS; the held will be amplihed by shearing induced by differential 
rotation inside the accreting torus. Assuming that the angular speed difference between 
the front and back edges nearly equals u, a held of initial strength B ^ 10^^ G will be 
amplihed to 10^^ G {i.e., well below equipartition with diherential rotation, an upper limit 
which it is not necessary to attain) in ^ 10 diherential rotations; the amplihcation timescale 
is ta ~ 0.01 s, comparable to the bursts’ rise-times. At this point the torus will radiate 
at the rate E ^ —B'^u^r^/Qc^ = 10^° erg enough to power a burst. The details of 
the conversion of Poynting hux into a magnetized relativistic wind have been discussed in 
both diherent and similar contexts by Usov (1992, 1994), Meszaros and Rees (1997b), Katz 
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(1997). Alternatively, one may invoke the Blandford-Znajek (1977) mechanism to extract 
the much larger hole’s rotational energy (Meszaros, Rees and Wijers 1998b). 

The energy release rate computed above of course exceeds the Eddington luminosity of 
the minitorus by more than 10 orders of magnitude. The energy deposition occurs outside 
the minitorus, but it is of course quite likely that some mass is blown off the torus, for 
instance by photons that hit it and heat it. This baryon outflow will certainly contaminate 
the zones of the relativistic wind that are most distant from the spin axis, because we do 
not expect any efficient angular momentum exchange in the process. Still, this implies that 
the outflow may be non-isotropic, and moderately beamed. The situation that we envision 
is similar to that described by Meszaros and Rees (1997b), with a smoothly tapering beam. 


3. A clean environment 


The volume immediately surrounding the SMNS will be vacated of baryons by two 
effects. First, the SN explosion accompanying the SMNS’s formation will sweep away ISM 
baryons. We can gauge the importance of this effect by considering the fact that most 
of the mass around the Crab pulsar, the only plerionic remnant for which relevant data 
are available, is well-known to be located in the hlaments, which amount to no more than 
~ 5Mq (Fesen, Shull and Hurford 1997), within about 3 arcminutes. At a canonical distance 
of 2 kpc, this corresponds to a present average baryon number density n ~ 10 cm~^. 
Assuming uniform expansion in the period from tsd to the Crab’s current age, implies that, 
at tsd, the average baryon density was n ~ 10^ cm~^. Within a distance D = 10^^ cm, more 
than enough to guarantee millisecond variability (which requires D k, 10^^ cm, Rees and 
Meszaros 1994) this corresponds to a total baryon mass M pa 3 x 10“^ Mq, below 10 “^Mq, 
the upper limit to get a signihcantly relativistic expansion. We regard the above estimate 
as an upper limit because the Crab is well-known to have an unusually low expansion 
velocity (a full factor of 3 below other Type II SNRs, Fesen, Shull and Hurford 1997) and 
because it considers the mass in the hlaments as if it were uniformly distributed over the 
whole volume rather than concentrated in the observed thin shell. 


The second effect is due to the large energy release by the magnetic dipole rotation 
after the explosion. 


E = - 




= —3 X 10"^^ erg s 
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since a large fraction of this energy will be converted to photons, the ensuing luminosity in 
fact exceeds the Eddington luminosity for the ~ 3Mq NS by about 4 orders of magnitude. 
This is of course identical to the well-known plerion model for SN remnants (Ostriker 
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and Gunn 1971, Reynolds and Chevalier 1984), except that the total energy released 
10 ^^ erg) is much larger; for this reason we cannot directly apply these models to 
our case, except to notice that its effect will surely be of lowering the naive baryon mass 
estimate obtained above 3 x 


4. Observational implications 

We now discuss observational signatures of the occurrence of this model. It shares 
with all models for which GRBs are generated when young, massive stars die a possibly 
exotic death, all predictions that associate GRBs with star forming regions. In particular, 
GRBs’ afterglows should be localized within star forming galaxies, and their redshift 
distributions ought to match closely the history of star formation in the Universe (Madau 
et ai, 1996, Lilly et ai, 1996). Furthermore, the optical afterglow may be at times 
absorbed by dust present in the star-forming regions, and the X-ray one may show 
absorption by large equivalent column depths by neutral elements. In fact, in the tsd 
years elapsed since the hrst explosion, the SN shock front will have reached a distance 
Rs = 10^® cm (uexp/3 X 10® cm s~^)(tsd/^0 yr). The total column depth of baryons, at 
distance is W = 10 ^^ cm“^(Mej /10 M 0 )(Rs/lO^® cm)~‘^ with the ejecta mass; this 
column depth may perhaps be observable (Murakami et ai, 1997). Thus we may sometimes 
expect to see the ionization edges discussed by Meszaros and Rees (1998a,b) due to dense, 
stationary outlying material. 

We argued above that in its maxi-Grab phase the SMNS will vacate of baryons a 
sizeable cavity around itself. Outside this cavity, however, there will be an outwardly 
increasing density gradient, a well-known property of Sedov solutions. In fact, matter 
behind the SN shock is subsonic with respect to the outwardly propagating shock, so that 
it must be in approximate pressure equilibrium with previously shocked material. However, 
since the shock is decelerating, just-shocked material will be at lower temperature than 
material shocked long before and, in pressure equilibrium, this implies higher densities. The 
exact Sedov solutions yield for the innermost regions uism oc where <7 = 3 /(7 — 1) = 9/2 
(Landau and Lifshitz 1979). If the time between the SN explosion and the burst is rather 
short, it may happen that the collision between different shells, which gives rise to the burst 
proper (Paczyhski and Xu 1994), takes place in a region which the maxi-Grab effect has 
not managed to empty of baryons, yet. In this case, the afterglow will propagate into the 
Sedov solution, i.e., inside an outwardly increasing density gradient. In the afterglow, the 
time tnr of the transition to the non-relativistic expansion regime is reasonably sensitive to 
the density gradient; Meszaros, Rees and Wijers (1998a) hnd tnr ~ where t^ is the 



duration the burst in the 7 -ray band, T ^ 100 is the initial Lorenz factor of the ejecta, 
and O' = (8 + 2d)/{3 + d), for adiabatic expansion (notice that we dehned as —d what they 
call +d). For expansion in a wind, d = —2, and thus O' = 4, while for expansion into an 
adiabatic Sedov solution, d = 9/2, q = 2.2, so that the transition to the non-relativistic 
regime is anticipated by up to four orders of magnitude in time. For = 1 s,q = 2.2, 
tnr ~ 7 hours. After tnr, the luminosity steepens considerably, to or steeper (Wijers, 
Meszaros and Rees 1997). This provides thus a natural explanation for those cases where 
either no afterglow is seen (in the hard X-ray GRB 970111, Feroci et ai, 1998; in the 
optical GRB 970828, Groot et al, 1998a), or it is seen to decrease more steeply than the 
other afterglows (GRB 980326, Groot et a/., 1998b). 

The formation rate of SMNSs in an L*-galaxy, l/tsMNS, is of course expected to be 
of order of, but less than, the rate of formation for pulsars, in our Galaxy currently about 
l/tp, with tp Ri 100 yr. If the formation rate is around a factor of 100 lower than this, as 
suggested by rough considerations on the range of mass and angular momentum over which 
SMNSs form, there is ample room to account for all observed GRBs, including a moderate 
degree of beaming. It may prove interesting to search for such objects. The fraction of 
the time between each birth through which they are active in the maxi-Grab phase is 
tsd/tsMNS, SO that there will be one L* galaxy every tsMNs/tsd containing an active one. 
Given the luminosity density of the Universe (jo = 1-7 x 10 ®Lq Mpc~^ in the V-band) 
and the typical luminosity of an L*-galaxy like our own, (L* = 10^° Lq in the V-band, 
Binney and Tremaine 1987) we hnd that the closest such active maxi-Grab should lie at a 
distance R = {3LRsmns/^T^ jotsdY^^ = 52 Mpc where we used Eq. E and tsMNS = 10 ^ yr 
which is necessary if SMNSs must account for GRBs. At this distance, the flux of Eq. 
corresponds to ~ 10“^° erg pulsed with frequency u = io/2ti Ri 1000 — 2000 Hz. 

If the Grab pulsar is anything to go by, we expect an X-ray flux of ~ 10“^^ erg 
(only a fraction of which will be pulsed) from a pointlike source well offset from the center 
of the galaxy. A detection of the pulsations from such a source is beyond the capabilities 
of Rossi XTE, but well within reach of very large area X-ray telescopes such as XMM and, 
especially, Gonstellation. 

In summary, we presented a SupraNova model for gamma-ray bursts, he., a model 
with a two-step collapse to a black hole not induced by accretion. The advantages of 
this model lie in its ability to keep the environment baryon-free, both because of the SN 
explosion that goes with the hrst collapse, which cleans up the surroundings, and because of 
the clean, silent, collapse of the NS once centrifugal support weakens critically. Second, we 
discussed a rotating magnetic dipole radiating much energy through its large rotation rate, 
which injects into the NS’s surroundings a luminosity exceeding Eddington’s by four orders 
of magnitude, and which speeds up as a consequence of energy losses. This too contributes 


to keeping the NS’s environment baryon-clean. A rough estimate suggests that the total 
environment pollution by baryons may be in the range of ~ 10“^ Mq. 
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